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Abstract 

We investigate a heterostructure system with a 
spin-orbit coupled semiconductor sandwiched by 
an s-wave superconductor and a ferromagnetic in- 
sulator, which supports Majorana fermions (MFs) 
at the superconducting vortex cores. We pro- 
pose a scheme of transporting and braiding the 
MFs, which only requires application of point-like 
gate voltages in a system with nano-meter pat- 
terns. By solving the time-dependent Bogoliubov- 
de Gennes equation numerically, we monitor the 
time evolutions of MF wave-functions and show 
that the braiding of MFs with non-Abelian statis- 
tics can be achieved by adiabatic switching within 
several nano seconds. 

Introduction. Majorana fermions (MFs) are 
particles equivalent to their own anti -particles .— 
Great research effort has been devoted to search- 
ing for MFs in condensed matter systems in the 
last decade since the particles can be used for con- 
structing topological quantum computers. 2 ^ Sys- 
tems that attract most attentions include Pfaffian 
v = 5/2 fractional quantum Hall state, 5 chiral 
p-wave superconductors (SCs), 6 one-dimensional 
(ID) spinless SCs, 7 superfluidity of cold atoms, 8 ' 9 
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to name a few. Recently, two heterostructure sys- 
tems are suggested as possible MF hosts, namely 
the topological insulator in proximity to s-wave 
SC— and the spin-orbit coupled semiconductor 
(SM) sandwiched by s-wave SC and ferromagnetic 
insulator (FI) (SC/SM/FI).-^ At the same time, 
a ID semiconductor nanowire with spin-orbit cou- 
pling in proximity to s-wave SC was investigated 
both theoretically and experimentally for realizing 
MFs.~~— Lately, a promising signal of MFs has 
been captured in the device of InSb nanowire/5- 
wave SC. 20 In these hybrid systems, the p-wave 
pairing are superseded by the interplay between 
proximity-induced s-wave superconductivity and 
strong spin- orbit coupling. 

It was revealed that MFs appear inside vortex 
cores in a chiral jc-wave SC, 6 and that non-Abelian 
statistics can be achieved by exchanging positions 
of vortices hosting MFs. 21 However, it is diffi- 
cult to manipulate vortices in experiments, which 
may hinder the realization of this genius idea. To 
circumvent this problem, MFs at sample edges 
of topological SC have been considered. 22 Mak- 
ing use of the topological property, edge MFs can 
be braided with the desired non-Abelian statistics 
merely by tuning point-like gate voltages on links 
among topological SC samples. In order to make 
the edge MFs stable, one needs to embed the de- 
vice into a good insulator. The size of topologi- 
cal SCs should also be chosen carefully since the 
wave-functions of edge MFs become too dilute for 
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Figure 1 : Schematic device setup for braiding MFs 
at vortex cores. There are four holes in SM layer 
(yellow platform) with one superconducting vor- 
tex (blue cylinder) pinned right beneath each of 
them. The electrodes at high-voltage states (pink 
rectangular prisms) prohibit electron hoppings in 
the regions below them, and thus connect effec- 
tively the holes; the blue rectangular prism denotes 
an electrode at zero-voltage state. 

large samples, which makes edge MFs fragile due 
to excited states. 

In this work we demonstrate that the core MFs 
can be liberated from vortices, transported and 
braided by applications of local gate voltages. This 
scheme takes full advantages of SC/SM/FI het- 
erostructure in the way shown schematically in 
Figure CD four holes are punched in the SM layer 
and one vortex is induced and pinned right be- 
neath each hole in the common superconductor 
substrate; three electrodes are placed above the 
small regions between holes, and the ones at high 
voltage state (pink rectangular prisms in Figure Q3 
connect holes by killing electron hoppings in SM 
locally. The linear dimension of holes and separa- 
tion between holes should be in the regime of tens 
of nano-meters: 22 for too large holes, the energy 
gap between the zero-energy MFs and the low- 
est excited states becomes very small, which lim- 
its the operation temperature, while for too small 
holes especially at the part of keyhole tail MF 
wave-functions interact with each other, which de- 
stroys the MF ground state. 

The key observation is that the geometric topol- 
ogy of the SM layer can be controlled by local 
gate voltages, and that when even number of holes 



are connected, core MFs fuse into quasi-particle 
states with finite energies, while one core MF 
exists when odd number of holes are connected, 
since the unified perimeter includes even and odd 
number of vortices in the two cases respectively. 
Core MFs can then be liberated from and trans- 
ported among vortices with a sequence of turning 
on and off gate voltages at the electrodes. In or- 
der to keep the topological protection, one needs to 
switch gate voltages in an adiabatic way. To moni- 
tor the dynamics of MFs, we solve numerically the 
time-dependent Bogoliubov-de Gennes (TDBdG) 
equation with high precision. It is then found 
that, with reasonable materials parameters, adia- 
batic manipulations of MFs and thus non-Abelian 
statistics can be achieved within several nano sec- 
onds, which reveals that the present scheme has 
the application potential for topological qubits and 
topological quantum computation. 

Topological superconductivity. We start from 
the tight-binding Hamiltonian of a SM with 
Rashba-type spin-orbit coupling on a square lattice 
in proximity to a ferromagnetic insulator 22 

#0 = - £ kflcCfc - 11 £ c^Cia + £ V z (cl t c n - c^qj 



+ £ Wa 



l,j,<7 1,(7 

C (i+x)a S y C i<*' ~ C (i+y)a S x C '^' 



+ h.c.\, (1) 



where fy and t l a are the nearest-neighbor hopping 
rates of electrons with reserved and flipped spin 
directions; /i and V z are chemical potential and 
strength of Zeeman field respectively; c. creates 
one electron with spin o at lattice site i; s = 
(s x , s y , s z ) are the Pauli matrices. 

The proximity-induced superconductivity in SM 
is described by 

tf sc = £(A44 + fc.c.), (2) 

i 

where A, is the pairing potential at site i. The BdG 
equation is then given by 



H A 

A f -6yH*o y 



vi>(?) = EV(r) (3) 



with the Nambu spinor notation ^(r) = 
[u+(r),u^(r),v±(r), — v-j-(r)] , which defines the 
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quasi-particle operator Y* = f dr[£ a u a (r)ca(r) + 
Vo(r)c a (r)}. 



Chern number 
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Figure 2: (a) Dispersions of two bands close to 
Fermi level (left panel) and Berry curvature Vx^/ 
(right panel) in the first Brillouin zone, (b) Sys- 
tem of two holes and two vortices with holes iso- 
lated (top panel) and connected (bottom panel), (c) 
Energy spectrum of several low-energy excitations 
at vortex cores when the holes are isolated. Up- 
per inset: energy spectrum of excitations at the 
sample edge. Lower inset: distributions of zero- 
energy quasi-particles. Results are for A = 0.5?o, 
V z = 0.8% /i = —AtQ and t a = 0.9?o with a sample 
of 400 x 200 sites, (d) Four spinor components of 
zero-energy states at the right hole, with the length 
and azimuth angle of vectors denoting the ampli- 
tude and phase of spinor fields respectively. 

To reveal the topological nature of the sys- 
tem, we evaluate the Chern number (also known 
as the TKNN index).— For this purpose, we 
treat first the system in k-space on the basis 

T 



•kt> C-k| 



reads 



The Hamiltonian then 



s(k) + V z s z +R(k)-s 



iA(k)s y 
-e(k)-V z s z +R(k)-? 

(4) 

where e(k) = —2to(cosk x +cosk y ) — fi, R(k) = 
(—2t a smky : 2t a s'mk x ). Diagonalizing H(k), 
we obtain the eigenstates <J>(k) and band struc- 
ture E(k) in the first Brillouin zone T 2 = 
[—n/a, n/a) ® [— n/a, n/a), as shown in Figure 
Ufa). The band topology is characterized by the 
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c = — A-(Vx 
2ni Jj* 



(5) 



with Berry connection for all occupied bands 

= E, J ;£ n (k)<o(*«(k)|4 M *n(k)) (/* = X,J). 

Integrating the Berry curvature V x srf shown in 
Figure |2ta) over the Brillouin zone, 24 we find 
that c = — 1 when the Zeeman field satisfies 
[Gu+4f ) 2 +A 2 ] 1/2 < Vz < (> 2 + A 2 ) 1/2 with 
fl < —2to, which is attributed to the topologically 
nontrivial energy gap at {k x ,k y ) = (0,0) while 
those at (0, —n/d), (—n/a,0) and (— n/a, —n/a) 
remain trivial.-^' 2 - 5 The nonzero Chern number 
indicates that the system is topological and thus 
may support MFs. The above result is consis- 
tent with the approximation treatment around T 
point- 
Core Majorana fermions. We study a finite 
sample with two separated holes in the SM and 
two vortices pinned right beneath them (see top 
panel of Figure Etb)). The typical size of a square 
sample is 300nmx300nm, which is divided into 
200 x 200 square grids, corresponding to Hamil- 
tonian matrix of dimension 10 5 x 10 5 in ©. By 
solving the BdG equation for this case, we obtain 
the energy spectra of excitations both at the vortex 
cores and the sample edge. Two zero-energy states 
are found at the holes, whereas no such state at the 
edge, as shown in Figure |2jc). We examine the 
four spinor components of the zero-energy states, 
and find u-\ = vl and = vT (displayed explicitly 
in Figure |2fd) only for the right hole), which re- 
sults in y> = 7, indicating that the two zero-energy 
states are Majorana states. 

It should be noticed that the excitation energy 
gap at vortices is about four times larger than 
that at the edge (see He)), which makes the core 
MFs more stable than edge MFs. 22 On the other 
hand, because the minigap associated with An- 
dreev bound states at superconducting vortex core 
proximity-induced in SM is roughly A 2 / (A 2 + 
ju 2 ) 1 / 2 ~ A with a small Fermi energy /! ~ A, 26 the 
influence from Andreev bound states to the core 
MFs can be neglected. It is in contrast to the case 
of SC exposed to vacuum where /! ^> A and thus 
the minigap is small in order of A 2 //i. 
Next, we impose a point-like gate voltage on 
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Figure 3: (a) Schematic diagram for braiding the two MFs at the left and right holes, colored by red and 
green respectively, (b) Distributions of the wave-functions of MFs obtained by solving TDBdG at several 
typical time steps. 



the region between the two holes to prohibit di- 
rect hopping of electrons by lifting the on-site en- 
ergy there as in the bottom panel of Figure 0b) 
(see also Figure [D, which merges effectively the 
two isolated holes into a unified one. Solving the 
BdG equation for this case, there is no zero-energy 
quasi-particle excitation state, since the combined 
hole includes two vortices.— 

Based on the above result, we can design a way 
to liberate and transport a MF from one vortex 
to another (see the first three columns in Figure 
|3]). Initially, the top and middle holes are con- 
nected together while the leftmost one is isolated 
and hosts a MF. We then combine these three holes 
by applying gate voltages on the region between 
the left and middle ones, which causes the MF to 
spread itself over the unified hole including three 
vortices (t = T). Finally, the MF is moved totally 
to the top by disconnecting the top hole from oth- 
ers (t = IT). It is noticed that the collapsing of 
MF wave-function on the top hole is a topologi- 
cal property, and is impossible for electrons. Solv- 
ing BdG equation numerically, we confirm that the 
gap between the zero-energy MF states and the ex- 
citations remains finite during the whole process, 
which guarantees the topological protection. 

Being able to transport one MF from one hole to 
another, we proceed to exchange positions of two 
core MFs in the system shown in Figure |3] (see 
also Figure [I])- Following the above transporting 
procedures, we further move the green MF from 
the right hole to the left one during t = IT ~ AT . 
At last, the red MF stored temporarily at the top 
hole is transported to the right one in the period 
t = AT ~ 6T . After the sequence of switching pro- 
cesses, the system comes back to the original state 



with the red and green MFs exchanged. 

In order to investigate the impact of the posi- 
tion exchange to the MF state, we monitor the time 
evolution of MF wave functions |*P(f)) by solving 
the TDBdG equation numerically 

ihj t \y(t))=H(t)\V(t)), (6) 

where H(t) is given in Figure |3]and depends on 
time via the hopping rates and t l a in the regions 
between holes, which are tuned adiabatically by 
the local gate voltages.— 

Even with the powerful computation resources 
available in these days, it is still hopeless to tackle 
this problem by directly diagonalizing the Hamil- 
tonian H(t) of dimension 10 6 x 10 6 for each time 
instant. Fortunately, it has been revealed that 
when the exponential operator is expanded by the 
Chebyshev polynomial 

exp[-iH(t)8t/h] = Y,c n (8t)T n (H), (7) 

n 

the coefficient c n (8t) decreases with n exponen- 
tially fast for small St.— Therefore, only sev- 
eral leading terms c n are necessary for a suf- 
ficiently accurate estimate of exp[— iH(t)8t/h]. 
Moreover, the Chebyshev polynomials can be con- 
structed based on the recursive relation T n (H) = 
2HT n -\(H) — T n -2(H), which reduces the com- 
putation time further. In this way, the time- 
dependent wave-function of the MFs can be ob- 
tained efficiently in an iterative way ^(t + 8t)) ~ 
exp[-iH(t)8t/h) 

The wave function = 0)) is the superposi- 
tion of the zero-energy states | ¥/,(())) and 1^(0)) 
at left and right holes respectively as shown in Fig- 
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ure Ha). We evaluate the projections of 
onto the initial states L = (*P/,(0) and 
0* = (*P#(0) and display them in Figure |4] 

for < t < 6T. The two MFs pick up opposite 
signs after exchanging their positions, which can 
be summarized by 

7l -> Yr, 7r -> -U- (8) 

Since ([8]) can be written as y L / R — >■ U^y L / R U with 
the unitary matrix [/ = exp(^y^7 L /4), the braiding 
of MFs satisfies non-Abelian statistics. 18,19,22,29 
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Figure 4: Projections of the MF wave-function 
\^{t)) obtained by TDBdG onto the initial states 
L = OP L (0)|¥(0) and R = (W R (0)\W(t)). 

The braiding rule d8) can also be understood by 
the three pairs of opposite motions of MFs across 
links between holes during the whole process dis- 
played in Figure |3l since each pair of opposite mo- 
tions contributes -1 and ( — l) 3 = — 1, which indi- 
cating that either Jl or Yr has to change sign. 22 

The above scheme for the braiding of MFs 
shows good scalability. For a device in Figure 
[5] constructed by the unit displayed in Figure [H 
one can exchange any two MFs by a sequence 
of switchings of point- like gate voltages. It is 
straightforward to figure out that any MF-pair ex- 
changing accompanies odd number of opposite 
motions of MFs across links between holes, which 
guarantees the non-Abelian statistics.™ Therefore, 
the present scheme has the potential to be used for 
large scale manipulation of topological qubits. 

Conclusion. We have shown that the Majo- 
rana fermions hosted by vortex cores in topolog- 
ical superconductors can be liberated from pinned 




Figure 5: Device array for the present scheme for 
manipulation of MFs. 

vortices, transported and braided over the pre- 
pared holes, taking advantages of the heterostruc- 
ture of 5-wave superconductor and spin-orbital 
coupling semiconductor. By solving the time- 
dependent Bogoliubov-de Gennes equation nu- 
merically, we monitor the time evolutions of Ma- 
jorana fermion wave-functions and demonstrate 
the non-Abelian statistics of adiabatical braidings 
of Majorana fermions. The present scheme only 
requires local applications of gate voltages, and 
minimizes possible disturbances to the Majorana 
fermions, which might be a challenging issue in 
other proposals based on end Majorana fermions 
in one-dimensional superconductors where gate 
voltages are necessary along the whole system. 
As compared with the edge Majorana fermions in 
two-dimensional topological superconductors, the 
core Majorana fermions are protected by a larger 
energy gap, and thus reduce the limitation on oper- 
ating temperatures. Therefore, the present scheme 
provides a more feasible way for manipulating 
Majorana fermions in large scale. 
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